def geometric_sum(N):
    a = 1/N
    r = 1/N
    s = 0
    for i in range(1, 1001):
        term = a * (r ** (i - 1))
        s += term
        if i % 100 == 0:
            print(f"The sum of the first {i} terms is {s}")
    return s


# Example usage
geometric_sum(2)
